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We study a unified mechanism for spin-density-wave (SDW) and superconductivity in a minimal 
model, in which itinerant electrons and local moments coexist as previously proposed for the iron 
pnictides [EPL, 88, 17010 (2009)]. The phase diagram obtained at the mean field level is in qual- 
itative agreement with the experiment, which shows how the magnetic and superconducting (SC) 
instabilities are driven by the critial coupling between the itinerant/localized electrons. The spin 
and charge response functions at the random phase approximation (RPA) level further characterize 
the dynamical evolution of the system. In particular, the dynamic spin susceptibility displays a 
Goldstone mode in the SDW phase, which evolves into a gapped resonance-like mode in the super- 
conducting phase. The latter persists all the way into the normal state above T c , where a strong 
scattering between the itinerant electrons and local moments is restored, as an essential feature of 
the model. 



I. INTRODUCTION 



High-temperature superconductivity found 1 -^ in the 
iron pnictides has attracted an intensive attention in re- 
cent years' The SC pairing of the electrons in these 
materials is less likely to be mediated by phonons, 
as suggested by the local-density-approximation (LDA) 
calculations 4 - as well as a variety of the experiments. 3 
The proximity of the SC state to the SDW phased in 
the phase diagram implies that the interplay between the 
magnetism and superconductivity might play an impor- 
tant role in understanding the pairing mechanism and 
other physical properties of the iron-based superconduc- 
tors. Combined with the high SC transition temperature, 
one finds an intriguing resemblance between this family 
of materials and the cuprates, in which superconductivity 
is generally believed of electronic origin. 

Nevertheless the electrons in the iron pnictides arc 
more itinerant than those in the cuprates, especially in 
the magnetic phase, where the electrons in the latter are 
localized due to a Mott transition to form local moments 
which become antiferromagnetically (AF) ordered at low 
temperature. By contrast, in the former, many differ- 
ent experiments have clearly demonstrated the itinerancy 
of the electrons in the SDW phase, including the multi- 
ple Fermi pockets as revealed by the angle-resolved pho- 
toemission spectroscopy (ARPES) measurements^— 
which are consistent with the LDA calculations jiir— the 
quantum-oscillation^ the transpor t 1 ' 2 ' 17 as well as the 
optical measurement^ Based on these experiments, one 
may reasonably view the magnetic order in this system 
as an SDW order formed by the itinerant electrons via 
Fermi-surface nesting] 13 ' 15 ' 19 From this point of view, it 
is natural to conjecture that the superconducting pairing 
is mediated by the collective magnetic fluctuations of the 
itinerant electrons. A lot of theoretical efforts have been 
done along this line, including the weak-coupling RPA 
theory ; 11 ' 12 ' 20 ' 21 the fluctuation-exchange approximation 



(FLEX) ; 22 ' 23 renormalization group (RG)j2£, functional 
renormalization group (FRG)^ 5 . - — and strong-coupling 
variational Monte-Carlo (VMC) 28 approaches. The sign 
change in the gap function between the electron and hole 
pockets is predicte d 11 ' 12 ' 25 as due to such unconventional 
magnetic origin of superconductivity. 

However, the magnetism in the iron pnictides has been 
also looked upon from the strong coupling side, where 
the electrons are localized via a multiband Mott transi- 
tion forming local moments^ - — Indeed, the magnetic 
ordering phase observed in experiment can be also de- 
scribed by utilizing a Heisenberg type J\-Ji modeh 34 ' 35 
A study of the SC state based on a doped Mott-insulator 
described by the t-J\-Ji model also yields 3 ^ a consistent 
pairing symmetry as compared to the experiment. This 
point of view is further supported by the first-principle 
LDA calculations which generally show the tendency for 
a large magnetic moment formation at low doping^ 7 - - — 
A strong experimental support for the local moment pic- 
ture comes from the observation of a linear temperature 
dependent magnetic susceptibility in a broad tempera- 
ture regime above the Neel temperature^ 2 - - — A recent 
neutron-scattering measurement 45 in FeTeo.35Seo.65 fur- 
ther observed a substantial magnetic moment persisting 
up to 300K above T c in the SC regime. 

To reconcile the two aspects of the itinerant and local- 
ized electrons exhibited in the iron-based superconduc- 
tors, a minimal model was proposed in Ref4&, in which 
the itinerant electrons and the local moments are con- 
jectured to coexist, based on the multiband nature of 
the system. Here the two separated degrees of freedom 
may be attributed^ 7 - to different 3d orbitals of the iron 
atoms, with the local moments formed via an orbital- 
selective Mott transition^ Recent dynamic mean-field 
theory (DMFT) calculations^ 9 - - — lend numerical sup- 
port for the possible orbital-selective Mott transition in 
the iron-based superconductors. Similar local-itinerant 
hybrid models with incorporating the detailed orbital 
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character s 52 ' 53 have been used to account for the mag- 
netic excitations observed in the iron pnictides. 

The most essential feature in the model of coexistent 
itinerant/localized electrons lies in the momentum match 
between the two degrees of freedom^ 6 - Namely, the char- 
acteristic wave vector of the magnetic correlation of the 
local moments, Q s = (7r, 0) or (0,tt), is commensurate 
with the typical momentum transfer between the hole 
and electron pockets of the itinerant electrons. Conse- 
quently the scattering between the two degrees of free- 
dom, which are coupled by the Hund's rule interaction, 
can get much enhanced in the normal state, rendering the 
system intrinsically unstable towards either the magnetic 
or superconducting ordering at low temperatures. At low 
doping, the magnetically ordered state can be obtained 46 
as composed of an SDW order of the itinerant electrons 
simultaneously locking with the collinear magnetic or- 
der of the local moments at the same wave vector Q s . 
Here the nesting effect in the itinerant degrees of free- 
dom alone or a pure J1-J2 superexchange interaction in 
the local moment part can be much weaker themselves 
in driving the magnetic transition. For example, a per- 
fect nesting of the Fermi surfaces can be easily removed 
by adjusting the chemical potential^ 6 - or by introducing 
a more realist band structure.— The pure collinear anti- 
fcrromagnetic (AF) ordering for the local moments can 
be also switched off by setting 4 - 6 J 1 — and/or making 
the subsystem in the disordered regime. But the hy- 
brid system can nevertheless experience an SDW phase 
transition due to the magnetic instability driven by the 
above-mentioned critical coupling between the two sub- 
systems. The residual scattering between the itinerant 
and localized electrons will become much reduced in the 
SDW state, which exhibit o 46 ' 52 a series of magnetic and 
charge properties qualitatively consistent with the iron- 
based superconductors at low doping. 

However, the SDW instability is only one of the pos- 
sible infrared fixed-points of this minimal model at low 
doping. With the increase of doping, the itinerant elec- 
trons can also form the Cooper pairs to gain the interact- 
ing energy between the two subsystems. By doing so, the 
superconducting energy gaps open up at the Fermi sur- 
faces of the hole and electron pockets, which also stabilize 
the system by reducing the residual scattering between 
the itinerant and localized electrons. Generally the ori- 
gin of the SDW and SC states can be explored on equal 
footing in the hybrid model of itinerant/localized elec- 
trons. 

In this paper, a systematic evolution of the SDW and 
SC states with doping is studied based on a simplified 
two-component model similar to the one proposed in 
Ref4£. Both the magnetic and SC instabilities are found 
in this model at low temperatures, and a phase diagram is 
determined at the mean-field level. In the SDW state, the 
dynamic spin susceptibility at the RPA level shows that 
the magnetic excitations are split into two branches com- 
posed of a low-lying Goldstone mode and a gapped high- 
energy mode dominantly contributed by the local mo- 



ments. The latter generally gets severely broadened due 
to the strong scattering between the itinerant electrons 
and local moments. On the other hand, the coherent 
Goldstone mode in the magnetic phase will be replaced 
by a low-lying gapped "resonance-like" mode when the 
system enters into the SC state, which can persist all 
the way to the high-temperature normal state, consistent 
with the neutron-scattering observations^ The band 
renormalization is also studied within the same frame- 
work. The results in this model study illustrate a phe- 
nomenology consistent with the iron-based superconduc- 
tors, in which neither the Fermi surface nesting for the 
itinerant electrons nor the superexchange interaction for 
the local moments play the direct role alone. Rather 
these effects get strongly enhanced via the Hund's rule 
coupling between the two subsystems with the momen- 
tum match. The latter effect makes the system generi- 
cally unstable against the SDW/SC orderings at different 
doping and therefore provides a unified mechanism to un- 
derstand both orders appearing in the iron pnictides. 

The remaining part of the paper is organized as fol- 
lows. In section II, we introduce the basic model, which 
includes a two-pocket description of the itinerant elec- 
trons and a nonlinear cr-model description of the local 
moments, which are coupled together by the Hund's rule 
coupling. In section III, we present a mean-field calcu- 
lation, which gives rise to a global phase diagram with 
both the SDW and SC orders identified at different dop- 
ings. Section IV is on the dynamic fluctuation beyond 
the mean-field approximation, including the RPA calcu- 
lations of the dynamic spin susceptibility, the uniform 
magnetic susceptibility in different phases and the band 
renormalization effect. Finally, the conclusion and dis- 
cussion are presented in section V. 



II. MODEL DESCRIPTION 

Our starting point is an itinerant-electron and local- 
moment hybrid model, previously proposed 4 ^ to describe 
the iron-based superconductors. We shall explore the 
emergent magnetism and superconductivity in this highly 
simplified model, which is to be specified below. 



A. Model Action 

The effective action includes an itinerant electron sec- 
tor Sit, a local moment sector Si oc , and a Hund's rule 
coupling term Sh as follows 

S c s = Sit + Sioc + Sh- (1) 
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1. Itinerant electron 

We consider a simple two-pocket model for the itiner- 
ant electrons, whose action reads 

St = X)4Hw + (2) 

k 

where k — (iu>, k) is the (fermionic) momentum- 
frequency vector, and = (cL,, Cpj_, c' M ^ t c' M ^) con- 
tains the creation operators of the itinerant electron 
for both spins (t and J.) and in the hole and elec- 
tron pockets denoted by T and M, respectively (see be- 
low). The momentum- frequency summation stands for 
J2k = Z 3 " 1 J2iu> Sfc> where /3 _1 = k B T, and this conven- 
tion will be adopted throughout this work. 

£fc is a 4 x 4 matrix which determines the band disper- 
sion. It may be written as 

£,k = -VPo - e kP3, (3) 

where pi = Oi ® cro is the 4x4 matrix defined by the 
Kronecker product of the Pauli matrices Oi (where ctq is 
the identity 2x2 matrix), [i is the chemical potential, and 
(k = k 2 / (2m) — €q models a parabolic band dispersion 
with the effective mass m and energy shift £q. 

According to ARPES, 6 ~ 9 magneto-oscillation 
experiments,— and LDA calculations,— ~— both Fermi 
pockets are small (with kp ~ n/10) and shallow 
(Ep ~ O.leV). So by setting the parameters at 
m = O.leV^ 1 and e = O.leV, we can produce a 
reasonable band structure as shown in Fig. [TJ 



band structure we focus on the electron pocket around 
(tt, 0) point and the hole pocket around (0, 0), which are 
connected by Q s = (ir,0). The scattering between (0,0) 
and (0, tt) pockets as connected by Q s — (0,tt) can be 
similarly treated. 

2. Local moment 

Consider an AF superexchange coupling Ji bridged by 
the As ions between the diagonal Fe sites. The local 
moments at the Fe sites may be divided into two sets of 
sublattices, each is described by an AF Heisenberg model, 
respectively. The effective Hamiltonian may be written 
as 

#ioc = ^ 2 ]T ]T M r -M r ,, (4) 

X=A,B {rr')eX 

where X = A,B labels the sites on different sublattices 
as shown in Fig. [2] and (rr') denotes the nearest neigh- 
boring sites in the same sublattice (or equivalently the 
next nearest neighboring sites in the original lattice). 
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FIG. 1: (Color online.) (a) The Fermi surface consists of a 
small hole pocket around the V point and a small electron 
pocket around M. Both pockets locate around the origin in 
the SDW Brillouin zone (SDW BZ). (b) Band structure. The 
nesting vector Q s = (tt, 0) connects the two bands. 

There are several aspects that we wish to comment on 
this simple two-pocket model. Firstly, at /x = 0, the hole 
and the electron Fermi pockets perfectly match under a 
momentum translation Q s = (tt, 0). In real materials the 
nesting is not perfect even in the undoped case,— _ — 1^— 
which may represented by a small but finite here. Sec- 
ondly, the detailed orbital characters are neglected with 
the pockets taken to be rotationally symmetric. Such a 
simple model can not account for some anisotropic phe- 
nomenon such as the nodal SDW gap as discussed in 
Ref^. Thirdly, for simplicity in the present two-pocket 



FIG. 2: Each circle in the figure represents a Fe atom. The 
Fe lattice is divided in to A and B sublattices. In each set 
of the lattice, the local moments interact by the nearest AF 
coupling (next nearest in the original lattice). 

The local moment M r will intrinsically fluctuate 
around the characteristic momentum Q s = (tt. 0) or 
(0, tt). So we may take M r = Mn r e t( ^ s ' r with n r a unit 
three-component real vector field, such that n r fluctu- 
ates smoothly in the space, which will be convenient for 
further field theoretical treatments. 

The low-energy AF fluctuation of the local moment in 
each sublattice may be described by a nonlinear er-model. 
The corresponding action reads Si oc = J2x=A b Sx, with 



Sx = T nx >-o (" 2 + ° 2q2 + nx " 



(5) 



where q — (iv, q) denotes the (bosonic) 
momentum-frequency vector, and nx.q = 
Jo ^ r St-gx e ^ q r ~ VT ^n r (T) is Fourier-transformed 
from the field n r in the X sublattice. rj 2 is the La- 
grangian multiplier that enforces the unit condition 
n 2 , = 1. According to the above definition, a long- 
wavelength limit at q — » actually corresponds to the 
real momentum — > Q s . 
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The coupling constant g and spin wave velocity c are 
related to the Heisenberg J 2 -model by g — 8J 2 , c = 
4J2M (with the Fe-Fe distance taken as the unit).— From 
the neutron-scattering experiments^ the typical spin 
wave velocity is around 0.3eV, so that we set c = 0.3eV, 
which in turn determines J 2 = 0.093eV and g — 0.75eV, 
assuming the magnetic moment M — 0.8 per Fe atom. 

By introducing a parallel field n = (jia + ub)/2 and 
an antiparallel field fi — (jia — tib)/^, the action can be 
further written as 



•Sloe =7T ^2 n_ q {v 2 + c 2 q 2 + i] 2 )r 
9 „ 



Yg^- q{V 



(6) 



3. Coupling term 

The itinerant electron spin Sy — S claCr at site r can 
be coupled to the local moment M r at the same site by 
a Hund's rule interaction 



-Jh ^ ' Mr ' S r) 



(7) 



where Jh denotes the strength of the effective Hund's 
rule coupling. Fourier-transforming to the momentum- 
frequency space yields the following action 



S c 



(8) 



fe . q 



where J = 2JhMS, M is the magnitude of the local 
moment, and S — 1/2 for the itinerant electron. In Eq. 
([5]), s = (si,s 2 ,s 3 ) with each Si (i = 1, 2, 3) as a 4 x 4 
matrix defined by s j = o\ ® Ui . 

The itinerant electrons will couple to the parallel field n 
which mainly causes the scattering between the T pocket 
and the M pocket at (tt,0), and to the antiparallel field 
h which mainly causes the scattering between the T and 
(0, 7r) pockets. 



B. Mass-gap equation 



where 770 is the solution of 77 at zero temperature, which 
reads 



2ttc 2 1 



3ff c 27 

where 7 = (g c /g)(l — Uq). The parameter g c — 27rc/(3A) 
is introduced to replace the momentum cutoff A which 
is needed to control the convergence of the momentum 
summation. 

The physical meaning of 77 is the mass gap of the 
bare spin wave excitation in the local moment sector. 
If g > g c , the mass gap will remain finite at zero temper- 
ature, indicating a disordered state with only short-range 
magnetic ordering for the pure local moment degrees of 
freedom. We shall see that due to the strong coupling be- 
tween the itinerant electrons and local moments, a true 
SDW order can be still induced, even if the bare local 
moments are in a disordered regime and do not order at 
T = by themselves. 



C. Propagators 



The itinerant electron single-particle propagator is de- 
noted by G(k) = — {ckc' k ), and represented by the ar- 
rowed line ■; in the Feynman diagram. G(k) is a 
4x4 matrix as c! has 4 components (2 pockets x 2 
spins). The propagator for the local moment is denoted 
by D(q) = — (n- q n q ), and represented by the dashed line 
===1 in the Feynman diagram. D (q) is a 3 x 3 matrix as 
n q has 3 components. 

The double lines in the Feynman diagram represent 
the dressed propagators, while the single lines, such as 

— ^— and , denote the bare propagators, Go and Do, 

respectively. The latter can be read out directly from the 
actions, S\t and S\ oc , as 



Go(fc) = ^-tkV 1 = 



and 



D ,u'(q) = 



(iu + fi)p Q - e k p 3 
(iu} + p) 2 -e 2 k 



> 2 + c 2 q 2 + r/ 2 



(12) 



(13) 



To simplify the calculation, we replace the Lagrangian 
multiplier r] 2 by its saddle point value. Evaluating the 
saddle point equation SS\ oc /S(rf) = at the one-loop 
level yields 



— {l-nl 
2<? V ° 



1 



v\ + c 2 q 2 + rj 2 



= 0, 



(9) 



where no = (ji q —o) denotes the mean- field value of the 
local moment. 

Carrying out the momentum and frequency summa- 
tions, rj can be solved from the above equation as follows 



"n = Vo + -5 In g i 1 + V 7 ! + -le-Pioj 



(10) 



III. MEAN-FIELD PHASE DIAGRAM 

In the effective action Eq. ([1]), even though the two 
subsystems, i.e., the local moments and itinerant elec- 
trons, may not be in a magnetic or SC ground state sepa- 
rately, the whole system coupled together will experience 
intrinsic magnetic and SC instabilities, which are studied 
on equal footing below at the mean-field level. 



A. SDW phase 

The SDW mean-field equation can be deduced from the 
following Dyson equation for the self-consistent Hartree 
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approximation 



(14) 



The SC mean- field equations are equivalent to 
the following Dyson equations of self-consistent Fock 
approximation 55 , 



G(fc) = G (fc) + G (fc)S d G(fc), (15) 
where the Hartree self-energy is given by 

S d = J$D 0>ii ,(q = 0) Sl £ Tr [G (k') a v ] . (16) 

k' 

On the other hand, the Hartree energy is related to the 
local moment mean field by 



— —JonoiSi 



(17) 



where n^i is the ith component of (ra g= o). This can be 

seen by separating out the q — term - Jo Et "o ' c k SCk 
from the summation in the coupling term Eq. ([8J . This 
term indicates that the local moment mean field no af- 
fects the itinerant electron Hamiltonian by adding the 
self-energy —JqUq ■ s, which should just be identified as 
the Hartree energy. The consistency between Eq. (fT^|) 
and Eq. (|17l) leads to the SDW mean field equation. To 
show this, we first take Eq. (1171) and evaluate the propa- 
gator from Eq. (ITSj) 



Ll \ K ) ~ 77. . i ,.\2 12 T2T2 ■ 



(iuj + fi) 2 — 
Substituting into Eq. (fH)|) yields 

4gJ$n 0i Si x - i 



72 ■ 

Jr. "■( 



o"-o 



E 



(ico + /1) 2 



"0 



After Matsubara frequency summation, 



S rf = — 



where El 



E 



2E k (cosh + cosh j3E k ) ' 



(19) 



(20) 



JqtIq. Comparing with Eq. (fTF)) . we 



arrive at the SDW mean-field equation 
4 ff Jq sinn ^^ fc 



ry 2 ^ 2£*(cosh fift + cosh /3£*) 



= 1, 



(21) 



from which the SDW order parameter no can be deter- 
mined for the given chemical potential and temperature. 



(22) 



The above diagrams correspond to 

G(k) = G (fc) - G (k)Z p (k)F(-k), 



(23) 



-F(k) = G (fc)£ p (fc)G(-fc), 
where the pairing energy is 

£„(*) = - J 2 D o,iv (k - k') Sl F (k') (-4) . (24) 
k> 

Mediated by the local moment fluctuation, the effective 
interaction between the itinerant electrons is described by 
the following action 



Smt - E 



,t J 

'k+p^-k 1 



k')c. 



k'Ck'+p 



k,k ; ,p 



where the vertex function reads 



T(<?) = -T = -yl>0,ii'(9) s i ® s i 



(25) 



(26) 



The eigenvalue of the kernel function T(k — k') stands 
for the effective pairing energy of the Cooper pair, whose 
form factor is given by the corresponding eigenvector. 
The most negative eigenvalue (hence the strongest pair- 
ing attraction) is found in the spin-singlet intra-pocket 
pairing channel with s -wave symmetry. In fact, simply 
by diagonalizing the matrix Sj ® Sj, it is easy to show 
that the greatest eigenvalues belong to the spin-singlet 
parings with opposite sign between the electron and the 
hole pocket. Among the spin-singlet pairing channels, 
the inter-pocket pairing would lead to the pocket-singlet 
which requires the gap function to be of p-wave symmetry 
and is not able to fully gap the Fermi surface, and thus 
the intra-pocket s^-wave pairing remains most favorable. 

By introducing a 4 x 4 matrix <v = 03 ® (202), the 
s^-wave paring operator can be simply denoted as 



k^ c -k — ( c Irt c - 



•(4 



kr^-kri c kri c -kr-t) 

i.Mt C -kM± ~ C kMi C -kMt 



(27) 



Therefore one may assume the pairing energy to take the 
same form 



B. SC Phase 

To deal with the SC order, we introduce the abnormal 
propagator for the itinerant electrons: F(k) = — (cfcC_fc), 
represented by a line with two arrows heading in opposite 
directions ■; . 



E p (fc) = A(k)d, (28) 

with symmetric gap function A(fc), i.e. A(— k) = A(fc). 
Although the other pairing modes may also appear in 
Tip, they are all omitted to simplify the derivation. In 
general, an s ++ -wave pairing can be also induced from 
s ± -wave pairing if the two Fermi pockets are no longer 
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symmetric in size (i.e. with a finite chemical potential 
/i in our model). However, according to our calculation, 
the s ± -wave will be the dominant component persisting 
up to larger /i. 

Then from Eq. (|23]i . one can find the solution of G(k) 
and F(k), 



where 



Z(k) 



and 



G(k) = 



Z(k) Po - E(k)p 3 
Z{k) 2 -E(k) 2 ' 



(iuj — fi) 



2 ^2 



A(fc) 2 



, , (iuj - fx) + 2/i, 
e fc 

2 A 



(iu — /i) 



(iu; — jit) 



= f(k)A(k)d, 



(29) 



(30) 



(31) 




0.00 0.02 0.04 0.06 0.08 0.10 



FIG. 3: (Color online.) The phase diagram calculated accord- 
ing to the mean field equations. The solid lines denote the sec- 
ond order phase transition boundaries, and the dashed lines 
denote the first order phase transition boundaries. The red 
point is a critical point. The SDW and SC phase boundaries 
are calculated independently, and the overlap region does not 
necessarily imply the coexistence of two orders. 



where 



/(*) = \ E 



2 ^ ( ; .^2 



' M 2 -(e fe + <;/i) 2 -A(fc) 5 



(32) 



Here we have projected out the components other than 
s ± -wave in the solution of F(k) as noted above. Substi- 
tuting the above solutions into Eq. (|2~4"|) yields 

E p (fc) = 3J 2 J2 A),n (* " *') / (*') A (k') d. (33) 
k' 

Comparing with Eq. (|28p , we arrive at the SC mean-field 
equation 

A(fc) = 3 J 2 D ,ii (k - k') f (k') A (k') . (34) 



This equation can be solved by numerical approach. 

To proceed with analytic analysis, we omit the k- 
dependence of A(fc) and replace it by a constant A. We 
also approximate Do(k — k') by its average value {Do) at 
zero frequency around the Fermi surface, 



(D ) = -g 



1 



c 2 (k- k'f + v 2 1 k k , eFS Vi 



= -\. (35) 



1 /2 

One finds rj 2 s = rj (rj 2 + Ac 2 kp) with k 2 F = 2me . 
Then the SC mean-field equation becomes 



1 



3ffJ 2 



(36) 



which, after the Matsubara frequency summation, yields 



1 = 



. 2 



where E 2 k = (e^ + qfi)~ + A 2 , which determines the SC 
gap A. 



C. Phase Diagram 

By solving self-consistently the two mean-field equa- 
tions, Eqs. (HJ) and ([37|. the phase diagram of the SDW 
and SC phases can be determined as shown in Fig.|3J In 
the following, we specify the choice of the parameters in 
the model. 

It is noted that there are three particular points in the 
phase diagram that can be calculated analytically based 
on the mean-field equations. They are the SDW critical 
temperature Tg DW at n = 0, the SC critical temperature 
Tj? at /i = (if not consider SDW), and the chemical 
potential /Usdwi a t which the SDW order disappears at 
zero temperature. The formulae read 



T S ° DW = 1.13T e 



-I/CA^Vsdw) 



T° = 1.13T e 



-l/(N°V S c) 



and 



Msdwi — Toe 



= T^-VC^Vsdw) 



(38) 



(39) 



(40) 



where N F — m/(2ir) is the density of state at the Fermi 
energy n = 0, T = e (A 2 /(2me ) - l) 1 / 2 , V S dw = 
igJo/r] 2 and Vsc = SgJo/Vefi- To be overall comparable 
to the experiments, we take Igbw — 150K and T° ~ 40K 



and make use of Eqs. (|38|) and ([39)) to fix the model pa- 
rameters: g c = 0.53eV and Jo = 0.39eV. Consequently, 
with all the basic parameters in our model given, the 
phase diagram is determined numerically as a function 
of the chemical potential in Fig. [3] 

Figure [3] shows a critical point at the SDW phase 
boundary, located at fi* = 0.014cV and T* = 76K, where 
the second-order phase transition boundary splits into 
two first-order phase transition boundary lines. Hence 
at T — 0, there exist two critical chemical potentials for 
the SDW transition: /isDWi = O.OlleV and /^sd\V2 = 
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FIG. 4: The SDW order as a function of chemical potential 
at T = 40K, showing the hysteresis loop of the first ordered 
transition with two distince transition chemical potential. 



0.022eV, respectively, as shown in Fig. 01 A similar first- 
order transition phenomenon has been also reported^ in 
some other theoretical approach based on a pure itinerant 
model. 

Figure[3j shows that the SC critical temperature is not 
sensitive to the chemical potential, in contrast to the 
SDW state, because the pairing mechanism of the present 
model is not sensitive to the Fermi surface nesting con- 
dition. On the other hand, the s^-wave superconductiv- 
ity does require finite Fermi surface densities of states 
in both hole and electron pockets to support the intra- 
pocket pairing. That explains why T c eventually vanishes 
at /j, = O.leV when the Fermi level touches one of the 
band bottoms of the hole/electron pockets. It is noted 
that beyond fj, — O.leV, superconductivity of other types 
of pairing symmetry is still possible in the present model, 
which will involve the intra electron or hole pocket pair- 
ing and require some incommensurate AF fluctuations of 
the local moment away from Q s . 

Therefore, the high-temperature normal state can be 
regarded as an unstable fixed point state in the present 
model, in which the itinerant electrons scatter strongly 
with the local moments due to the momentum match at 
Q s and the model defines the relevant degrees of freedom 
that render the system flow into either an SDW or SC 
ordered phase, depending on doping, as the temperature 
lowers. The mean- field equations Eqs. ([21]) and (l37l) de- 
scribe, respectively, how the SDW and SC orders emerge 
from such a normal state, with the phase diagram in 
qualitative agreement with the iron superconductors^ 

Finally, we point out that although Fig. [3] suggests 
that the SDW and SC phases may coexist at low doping, 
in mapping out the phase diagram in the figure, only 
the maximal temperature of Tsdw & n d T c is shown at a 
given fi with assuming the vanishing of the other order in 
the mean- field equations. In other words, in order to de- 
termine the coexistent SC state inside the SDW regime, 
one needs to further incorporate the detailed competi- 
tion of the two orders into the self-consistent mean-field 
equations, which can be straightforwardly done by gener- 
alizing the above formulation. But this is not considered 
here not only for the sake of simplicity, but also because 
we wish to emphasize that the mutual interplay between 
the SDW and SC orders are not essential in driving their 



own formations in our model. 



IV. DYNAMIC FLUCTUATIONS 

The interaction between the itinerant and localized 
electrons has played a crucial role in resulting in the SDW 
and SC states, as described by the global diagram in Fig. 
[3] In the following we further investigate the evolution 
of dynamic fluctuations beyond the mean-field approxi- 
mation in these phases. 



A. Dynamic spin susceptibility 

To study the low-energy spin dynamics around the 
SDW wave vector Q s , we first consider the RPA cor- 
rection to the propagator D(q) of the local moment by 
the Dyson equation 



or 

D(q)=D (q)+D (q)U(q)D(q), 
which is solved formally as 

D(q) = D0{q) 

{q > i-D ( q )n{ q y 

where the RPA bubble II is given by 



(41) 



(42) 



(43) 



IV (ff)=-0= J " E Tr & G (kWG(k + q)} . (44) 



Here as the propagator of the itinerant electron, the par- 
ticular form of G(k) depends on the mean field states: 
in the normal state, it takes the form of Eq. (fT2j) ; in the 
SDW state, it is given by Eq. (18]). While in the SC state, 
the contribution from the abnormal propagator should be 
included as well 

n H '(q) = - 

=J^Tr[ Si G(k)s if G(k + q)] (45) 

k 
k 

where Eqs. (I2U1) and (f3"Tj) are used. 

By noting — jQ 2 H(q) represents the spin susceptibility 
of the itinerant electrons at the mean-field level, one can 
similarly write down the spin susceptibility of the itin- 
erant electrons at the RPA level, and finally obtain the 
following total spin susceptibility 



D (g) + J ' 2 n( g ) 
i-A)(e)n(<?) 



(46) 
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The inelastic neutron-scattering spectroscopy (INS) 
can measure the dynamic spin susceptibility as the imag- 
inary part of x{q)j ~ ^0-X( v + *0+, q), obtained after the 
Wick rotation iv — > u + i0+, which is presented in Fig. [5] 
in different phases (see the figure caption for the details) . 
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FIG. 5: (Color online.) The calculated spectral function of 
dynamic spin susceptibility. Darker shade indicates higher 
intensity, (a) in the normal phase (fi = 0.03eV, T = 50K). 
(b) in the SC phase (jj, = 0.03eV, T = 30K). (c) and (d) 
are both in the SDW phase (fi = O.OOeV, T = 80K). (c) 
shows the spectrum of transverse fluctuation, and (d) is for 
the longitudinal fluctuation. The red dashed curve marks out 
the bare spin wave dispersion. 



In the normal state, as Fig. 5(a) shows, the spectrum 
of the local moment fluctuation becomes very fuzzy when 
it immerses into the continuum of the itinerant electrons 
as indicated by the dome-shaped shadow area around 
q < 0.2. For comparison, the bare dispersion of the local 
moment spin wave is marked out by the dashed curve 
(which is gapped as in a disordered regime of the nonlin- 
ear a- model as noted before). The smearing of the spec- 
trum is clearly the result of the strong scattering between 
the itinerant and local moment in the region around Q s . 
At q ~ 0, a hot sp ot at the frequency slightly below 2\i 
can be seen in Fig. 5(a) This is because to create a spin 
flip at q = involves a pair of electron and hole excita- 
tions at the r and M pockets respectively which costs at 
least energy 2/i to go across the Fermi surface, and on 
the other hand, the gap of the local moment fluctuation 
is higher than this energy such that the scattering dimin- 
ishes. As a matter of fact, such a "resonance- like" mode 



becomes even sharper in the SC phase (Fig. 5(b) I simply 
due to the further reduction of the scattering with open- 
ing the SC gap. It may account for the "resonance mode" 
found in the INS experiment, 57 which indeed persists all 
the way to the normal state. 

In the SDW state at low doping, inside the SDW gap 
of the itinerant electrons, the fuzzy continuum is replaced 
by some emergent collective modes. The transverse spin 



fluctuations (in the directions perpendicular to the or- 
dering direction) become the gapless Goldstone modes 
(Fig. 5(c)), which is consistent with the previous RPA 



calculation using a more complicated five-band model for 
the itinerant electrons^ On the other hand, the longitu- 
dinal fluctuation (along the ordering direction) remains 
gapped as shown in Fig. |5(d)| 

The existence of the Goldstone mode can be proven 
rigorously at the RPA level. Since the RPA bubble in 
the SDW phase has a rather simple expression at zero 
frequency and momentum 



n^(o) = -4j 2 ^(^ - 



Jonoin i> 



sinh j3Ek 
2_Efc(cosh Pfi + cosh/3_Efc) ' 



(47) 



Let us suppose that the SDW ordering is along the 3rd 
direction in the spin space, i.e. n.o,i = ^o,2 = and 
n ,3 + 0. Then n n (0) = II 22 (0) < h 33 (0) < 0, and by 
comparing with Eq. (|20l) and referring to Eq. (|13l) , it is 
recognized that 



Srf = ^«n ll( o) 



-JonoiSiAuiWlMO). (48) 



As the self-energy is determined self-consistently from 
the SDW mean-field equation E^ = — JonojSj, at the 
mean-field saddle point, we have D o ,ii(0)IIii(0) = 1, 
which leads to a pole of D\\ at q = according to 
D = (1 — DoH)^ 1 Dq, proving the existence of a zero 
energy collective mode, i.e. the Goldstone mode. The 
same argument applies for the D22 component as well. 
Also taken into account the fact that D12 = 0, it can be 
concluded that there are two Goldstone modes, both are 
in the transverse directions. As for the D33 component, 
since H33 > IIu such that A),33(0)n 33 (0) < 1 (note that 
D o (0) is negative), no pole can appear at q = in gen- 
eral, meaning that the longitudinal mode is still gapped. 



B. Uniform Susceptibility 

The total uniform susceptibility x — Xioc +Xit at fi = 
is presented in Fig. 6(a)| in which the contributions from 
both the local moment and the itinerant electron degrees 
of free dom, i.e., xioc and Xitj are also given, respectively. 
In Fig. 6(b)] the uniform susceptibility is shown at dif- 
ferent /x's where the low-temperature phases are cither 
SDW or SC. 

To probe the uniform susceptibility for the local mo- 
ment, we add a Zeeman term — M J^- h ■ n% to the local 
moment Hamiltonian, where h is the uniform magnetic 
field. Then the local moment action is modified from 
Eq. ([6]) by the replacement v — > v + imh with m = 0, ±1 
denoting the quantum numbers of the three spin wave 
modes respectively. By integrating out the local moment 
degrees of freedom (i.e. n q and fi q fields), the free energy 
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for the local moment reads 



^2 ^^(v + imh) 2 + c 2 q 2 + i] 2 ). 

m=0,±l q 



(49) 



Then the uniform susceptibility can be obtained from the 
second order derivative xioc = —d 2 F\ oc /dh 2 taken in the 
h — >• limit, 



yi - n o i iV -v 2 + c 2 q 2 + v 2 

Aloe ~ / j t 1 t 2 2 i 2\? ' 

g ' (^ z + c^g^ + ify 



(50) 



Carrying out the frequency and momentum summation, 
we get 



1 



Xloc 



g tt(3c 



2 



(51) 



where the function Y(x) = xcothcc — ln(2sinha;). 

In the high temperature limit, according to Eq. (|10p . 
the spin wave mass gap rj increases linearly with temper- 
ature as r\ = 2ksT, then the function Y tends to a finite 
limit Y(l) — 0.458, resulting in a linear-T behavior 



k B T 

Xioc = 0.458 



(52) 



which will dominate the total uniform susceptibility at 
high temperature, consistent with the experiments . 42 ' 43 
On the other hand, the itinerant electron uniform sus- 
ceptibility can be evaluated from 



Xit =-S 2 Y / Tr[ Ti G(k)T V G(k)}, (53) 



k 



where S = 1/2 for the itinerant electrons, and the matrix 
Tj = (To <g> <7j represents the spin operator. The particular 
form of the propagator G(k) will depend on the order 
in the itinerant electron state. In general the frequency 
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FIG. 6: (Color online.) (a) The uniform susceptibility v.s. 
temperature at /i = 0. The dashed lines indicate the contri- 
butions from either the local moment (blue) or the itinerant 
electron (red). The solid line is the total uniform suscepti- 
bility, (b) The total uniform susceptibility curve at various 
chemical potentials [i, all shows a rapid drop in the ordered 
phase at low temperature. 



summation involved can be complicated. However, to the 
leading order of approximation, we have 

Xit =-l^ n 'F( E i)+ n 'F(E k )} + O(n 2 ,A 2 ), (54) 

k 

where n' F is the first order derivative of the Fermi dis- 
tribution function, and E k provides the band struc- 
ture. For normal state E k = ±e k — (a, for SDW state 



±(el + J 2 no) 1/2 ~ Mj and for sc state E t = 
((±efe — /i) 2 + A 2 ) 1 / 2 . The remaining terms are of the 
second order of the order parameters. Since the function 
n' F (E k ) peaks at E k = 0, so if the itinerant electron 
band is gapped from the Fermi surface, the contribution 
to the uniform susceptibility will decrease rapidly, which 
accounts for the quick drop the total uniform suscepti- 
bility in the ordered phase in Fig. |6(b)} 



C. The renormalization of the itinerant electron 
band 

Now we consider the self-energy correction due to the 
scattering of the itinerant electrons with local moments, 
which is given by 



E e (fc) = 



= - = J$J2 D n'(k' - k) 3l G{k') Sl ,. (55) 



Here the local moment propagator D is taken from Eq. 
(j43| as the RPA-corrected one, while the bare single- 
particle propagator G is given by Eq. (|12[) in the normal 
state and Eq. (fT8"|) in the SDW state. The renormalized 
single-particle propagator obtained from Dyson's equa- 
tion 



G(k) = [G(k)- 1 - S e (fc)]" 



(56) 



determines the spectral function after a Wick rotation to 
the real frequency domain by 



A(u, fc) = -2ImG(cj + i0+, fc). 



(57) 



The result for the hole pocket around the T point is shown 
in Fig. [3 The pocket is slightly more shallow in both the 
normal and SDW phases, compared to the bare disper- 
sion as indicated by the red dashed curve. 

This band renormalization effect can be understood by 
looking at the frequency dependence of the momentum- 
accumulated self-energy £ e (cj) = J2 k S e (cj, fc), as shown 
in Fig. [8] The negative imaginary part typically has a 
valley shape, due to the reduced scattering rate within 
the local moment gap ±7/. It can be well approximated 
by — 2ImE e (tj + i0+) cx uj 2 for small frequency uj. Ac- 
cording to the Kramers-Kronig relation, the real part 
of the self-energy should follow Rc£ e cx -co, meaning 
that the self-energy correction reduces the electron en- 
ergy above the Fermi level and increases it below the 
Fermi level, thus always squeeze the electron pockets. 
This partly account for the reduced pocket depth gener- 
ally observed in ARPES experiments^— compared to the 
LDA calculations 1 ^—. 
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FIG. 7: (Color online.) The calculated itinerant electron spec- 
trum for the r pocket in (a) the normal phase (fi = 0.03eV, 
T = 70K) and (b) the SDW phase (fx = O.OOeV, T = 70K). 
Darker shade indicates higher intensity. The red dashed curve 
marks out the bare dispersion of itinerant electron. 
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FIG. 8: (Color online.) The momentum- accumulated self- 
energy E e (ct)) in (a) the normal phase (fi — 0.03eV, T = 70K) 
and (b) the SDW phase (jj, = O.OOeV, T = 70K). The red solid 
curve represents the real part ReE e , while the blue dashed 
curve represents the negative imaginary part — 2ImE e . 

V. CONCLUSION 

In this paper, we have presented a systematic study 
of the itinerant electron and local moment hybrid 
model 46,52 for the iron-based superconductors. The mi- 
croscopic origin of both the itinerant electron and local 
moment degrees of freedom are all from the 3d-orbitals 
of the iron atoms. As a renormalization flow at low en- 
ergy, part of the 3d electrons is conjectured to form lo- 
cal moments through an orbital-selective Mott transition. 
Here as a simplification, a two-pocket band structure is 
adopted for the itinerant electrons without considering 
their Coulomb interaction. A robust short-ranged AF 
fluctuations around the momentum Q s is incorporated 
for the local moment part via a non-linear a description 
tuned in a disordered regime, which can persist up to 
high temperature in the normal state. Thus, in this min- 
imal model, an SDW/AF instability in either subsystem 
is not intrinsically present when they are decoupled. A 
Hund's rule ferromagnetic interaction then couples these 
two subsystems together. 

What we have established in this work is that such a 
simple model is generically infrared-unstable against ei- 
ther magnetic or SC ordering at low doping, thanks to the 
"resonant" scattering of the itinerant electrons between 
the hole-electron pockets by the local AF fluctuations 



of the local moments around Q s . In other words, the 
itinerant electrons form an SDW/SC order by a strong 
coupling to a background AF fluctuations of the pre- 
formed local moments with a momentum match. The 
phase diagram in Fig. 3 is qualitatively in agreement 
with the experimental ones, in which the Cooper pairing 
is not glued by the Fermi-surface-nesting driven collec- 
tive fluctuations of the itinerant electrons which would 
otherwise result in a much weaker pairing strength in a 
much narrower doping regime, close to the SDW phase, 
than what has been shown in Fig. 3. The effective glue 
provided by the magnetic fluctuations of the local mo- 
ment automatically favors the s^-wave paring symme- 
try here. The presence of the local moments further ex- 
plains the high-temperature linear- T dependence of the 
uniform magnetic susceptibility (Fig. 6) in the normal 
state. In particular, the strong scattering between the 
itinerant and localized electrons is represented by the 
dynamic spin susceptibility shown in Fig. 5, which il- 
lustrates how the Goldstone mode in the SDW state be- 
comes a "resonant-like" mode in the SC state as well as 
its evolution in the normal state. The strong signature 
of the itinerant /localized electron coexistent picture seen 
in Fig. 5, including both low and high energy parts, can 
serve a very useful qualitative prediction for the neutron- 
scattering measurement even if the comparison may not 
yet be quantitatively due to the highly simplified nature 
of the model. 

Therefore, the minimal model studied in this paper 
may be generally used to describe the low-energy physics 
in a multiband electron system in which the electrons 
in some more localized orbitals may first form short- 
ranged (fluctuating) SDW order at a higher character- 
istic temperature (called the hidden local SDW order in 
Ref. 33 ). Then at lower temperatures, the electrons in 
more itinerant orbitals can be naturally driven into a 
true SDW order or SC state via the Hund's coupling to 
such a preformed local SDW background. In contrast to 
the scenario^ that an electron may carry both a coher- 
ent itinerant and an incoherent local moment signatures, 
in analog to a single band case at an intermediate cou- 
pling, the multiband case provides with us an alternative, 
but simpler possibility, i.e., via the orbital-selective Mott 
transition, itinerant and localized electrons may be ex- 
plicitly separated as independent degrees of freedom. 
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